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Abstract. We discuss here various methodologies and an 
optimal strategy of the temperature and emission measure 
diagnostics based on Hinode X-Ray Telescope data. As an 
example of our results we present determination of the tem- 
perature distribution of the X-rays emitting plasma using fil- 
ters ratio method and three various methods of the calculation 
of the differential emission measure (DEM). We have found 
that all these methods give results similar to the two filters 
ratio method. Additionally, all methods of the DEM calcu- 
lation gave similar solutions. We can state that the major- 
ity of the pairs of the Hinode filters allows one to derive the 
temperature and emission measure in the isothermal plasma 
approximation using standard diagnostic based on two filters 
ratio method. In cases of strong flares one can also expect 
well conformity of the results obtained using a Withbroe - 
Sylwester, genetic algorithm and least-square methods of the 
DEM evaluation. 



1 Introduction 

The solar corona with temperatures greater than a million 
degrees emits mainly in X-ray domain of the spectrum. Its 
emission contains both the continuum from thermal plasma 
and line emission of highly ionized elements. The earliest 
phases of solar flares are dominated by non-thermal emis- 
sions (hard X-rays), the remainder of the events manifests 
itself primarily by thermal emissions at wavelengths ranging 
from X-rays through visual. The thermal radiation comes 
from plasma having a wide range of temperatures. 

Temperature and emission measure diagnostic of the ac- 
tive region and flaring loops is very important to model the 
loop structures, and analysis the heating sources and their 
(heating) distribution along the loop. 



For the optically thin emission, assuming isothermal 
plasma temperature, the well known method of filter ratio, 
based on ratio of the obs erved fluxes in two s elected energy 



(or wa velengths) bands (iVaiana et al.1 (1 19731) . iThomas et al 
(1985)) can be used. In the case of multi-temperature 
plasma the filter ratio method derives the line-of-sight av- 
erage isothermal temperature and emission measure. There- 
fore, in this approach it is more adequate to use many filters 
or line in terms of d ifferential em i ssion measure distribu- 
tions DEM (T) /dT dPottash et al.1 d 1964b . ISylwester etaL 
(1980)). Extended discission of the DEM versus isothermal 



approximation can by found e.g. in Sch melz et al. (2001) 
Cirtain et al.1 d2007l) and references therein. 
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Hinode is a new advanced optical, UV and X-ray solar 
satellite launched on 22 September 2006. It is equipped 
with three instruments: Solar Optical Telescope (SOT), 
EUV Imaging Spectrometer (EIS) and the X-Ray Telescope 
(XRT). XRT is a high resolution Wolter-I type grazing inci- 
dence telescope, which is sensible to the emission of plasma 
in the temperature range 5.5 < log T < 8. Such wide temper- 
ature coverage is realized by nine X-ray filters, each with its 
own passband and therefore a different response to plasma 
temperature. These filters are mounted on two filter wheels 
which allow many combinations of filters and thus can be 
used for a temperature diagnostic. Figure[T]presents the the- 
oretical temperature response functions of the ten XRT filters 
combination for unit volume emission measure of 10 44 cm~ 3 
for each filter used in this paper. XRT has 2048x2048 back- 
illuminated CCD array with 13.5 /im pixel size. It allows the 
Sun soft X-ray images to make in 34 x 34 arcmin FOV with 
spatial resolution of 1 arcsec. Exposure times are in range 
from 1 ms to 16 sec. Detailed description of th e instrument 
can be found in the paper by lGolub et al.1 (120071) . 

In this paper we will check if XRT response characteristics 
allow diagnostics of temperature (T) and emission measure 
(EM) of the coronal plasma. We also describe the obser- 
vations selected for the analysis. Additionally we present 
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Fig. 1. Theoretical temperature response functions of the selected 
XRT filters. 



and discuss filters ratio method in isothermal approximation 
and three different methods for differential emission measure 
(DEM) analysis in multitemperature approximation, respec- 
tively. 

In all calculations we used the temperature response func- 
tion as given by SolarSoft and presented in Figure Q] 

2 Data analysis 

For our work we selected three typical X-ray emitting ac- 
tive regions observed by Hinode/XRT with 6 or more dif- 
ferent filters. To prepare the selected data we used standard 
SolarSoft routine xrt.prep.pro which is similar in nature to 
sxtjprep.pro (Yohkoh software). This routine converts raw 
data into index and structure arrays. Additionally, this pro- 
cess removes several instrumental effects. We could men- 
tion some of the most important, for example: replacement 
of the near-saturated pixels for value greater than 2500 DN 
with 2500 DN, remove radiation-belt/cosmic-ray and streaks, 
calibration for read-out signals, removing of the CCD bias, 
calibration of the dark current, normalization of each images 
for exposure time. Using processed data we calculated their 
differential emission measure distributions and using two fil- 
ters ratio method, temperature and emission measure maps 
in isothermal approximation. 

2.1 Two filters ratio diagnostic 

The filter-ratio method is a technique commonly used for 
plasma diagnostics based on data from broad-band X-ray 
imaging systems. This method requires analysis of two or 
more intensity maps (images) of the investigated structures 
obtained with various broad-band filters. It delivers maps of 
temperatures and line-of-sight emission integrals or emission 
measure distributions as a function of temperature. 



For a multi-band instruments like XRT, under assumption 
that the emitting plasma is optically thin, the observed flux 
Ff. in filter k in a single pixel is: 



F k = 



f(T, X)r] k (X)EM(T)dXdT 



(1) 



TJX 



where f(T, A) is coronal plasma emissivity as function of 
the plasma temperature T and the wavelength A, calculated 
at distance of 1 AU using CHIANTI spectral code (including 
both emission lines and continuum). f] k (X) is an effective 
area of the instrument which includes the geometric area and 
the reflectance of the telescope mirror, the entrance and the 
focal plane filter transmissions and the CCD quantum effi- 
ciency. EM(T) is a collection of emission measures over a 
range of temperatures located along the line of sight. Defin- 
ing temperature response function Pk(T) of filter A: as a con- 
volution of the instrument effective area r\ k and emissivity of 
solar plasma /: 



Pk = J f(T,\)n k (X)dX (2) 
we can rewrite Equation 1 as: 

F k = J P k (T)EM(T)dT « £ P k {T t )EM{T t ) (3) 



For isothermal approximation this Equation can be written 



as: 



F k =P k (T)EM(T) (4) 
Similar to lThomas et al. ( 1985 ) using Equation 4 one can 



calculate ratio R of the observed fluxes (F) in two filters la- 
beled "1" and "2" for a co-aligned pixel area. This can be 
written as: 



R{T) ~T 2 -pW) 



(5) 



According to Equation 5, the ratio R of observed fluxes de- 
pend on plasma temperature T only. In the next step, one can 
calculate the emission measure EM(T) = F k /P k (T). 

This procedure can be used on all combinations of the 
XRT filters, but it has a limited usability. Especially, ratio R 
can be sometimes an ambiguous function in the whole range 
of temperatures. 

Using procedure described above we made temperature 
and emission measure maps for an area limited by an isoline 
0.3 of the brightest pixel on the image. 

2.2 Analysed events 

The first analysed region was a relatively faint active region 
NOAA 0955. Observations were made on 13 May 2007 at 
around 18:00 UT. Images taken in ten filters combination are 
available. Temperature and emission measure maps were cal- 
culated using two pairs of filters AL_poly / ALJhick and 
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Fig. 2. Temperature and emission measure maps calculated from two niters ratio, AL_poly/AL_thick for the active region NOAA 0955 
observed on 13 May 2007 at 17:57:30 UT. 
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Fig. 3. Temperature and emission measure maps calculated from two filters ratio, AL_poly/Be__thin for the active region NOAA 0923 
observed on West solar limb on 20 November 2006 at 21 : 16: 19 UT. 



Ti__poly / AL_thick. As an example the resulted maps for 
AL_poly / Al__thick pair are presented on Figure [2] For both 
pairs of filters the filter ratios are sensitive for temperatures 
from log(T) — 6.1 up to log{T) — 7.7. However, for all 
structures of this quiet active region we obtained tempera- 
tures in ranges 2.0-3.3 MK and 1.8-2.4 MK for the AL.poly 
/ AL_thick and Ti__poly / Al__thick pairs respectively. The 
range of the observed temperatures is typical for solar quiet 
active regions (see e.g. iGolub and Pasachoff ( 1997 )). 

An active region NOAA 0923 was seen on the edge of 
the West solar limb between 19 and 23 November 2006. 
Observations analyzed here were made on 20 November 
2006 at around 21:16 UT. Images taken in six filters are 



available. Temperature and emission measure maps were 
calculated for two pairs of filters: Al__poly/Al_mesh and 
Al__poly/Be__thin. The resulted maps for Al__poly/Be_thin 
filters pair are presented on Figure [3] The highest temper- 
atures range of 3.5 MK are observed in the brightest loop, 
seen just on the solar limb. 

The same active region NOAA 0923 was observed one 
hour later with seven filters. Temperature and emission 
measure maps were calculated for two pairs of filters: 
Al__poly/Al__mesh and Al__poly/Be__thick. The resulted 
maps for Al__poly/Be__thick pair are shown in Figure [4] 

Due to the construction of the filter wheels, XRT images 
in various filters are only obtained consecutively. Thus an 
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Fig. 4. Temperature and emission measure maps calculated from two filters ratio, AL_poly/Be__thick for the active region NOAA 0923 
observed on West solar limb on 20 November 2006 at 22:01:25 UT. 



important assumption underlying the filters ratio diagnostic 
(and also DEM analysis described below) is that there was 
no substantial change in the morphology and energetics of 
the event between two images using various filters. This as- 
sumption is usually fulfilled in the case of quiet Sun and ac- 
tive region observations. Checking the light-curves for dif- 
ferent filters we proved that this assumption is also true in 
the observations above described. Moreover, in the case of 
small B 1 .5 flare images with different filters were obtained 
first around the maximum and the second one on decay phase 
of the event. We use this pair because the second filter was 
only one "thick" filter available. The exact co-alignment of 
images is likewise important, especially for the DEM calcu- 
lations in the individual pixels. 

Using 'harder' ALpoly/Be_thick filters pair instead of 
'softer' Al_poly/Be_thin pair results in obtaining much higher 
temperatures. Temperatures in the brightest loop reached av- 
erage values up to 13 MK. Such a high temperature could be 
connected to the appearance of a small B1.5 class solar flare 
observed in this active region. The flare started at 22:00 UT, 
reached its maximum at 22:05 UT and ended at 22:09 UT. 
In fact temperature obtained from GOES data for this small 
flare reached the value around 10 MK. 

It is possible with XRT to register images of reasonable 
S/N with cadences 1 or 2 seconds. In the case of solar flares 
where fluxes can vary rapidly it is possible in principle to 
interpolate the intensities of one image to match the time at 
which the image in the other filter was acquired. Examples of 
temperature and DEM analysis of solar flare are well kno w 
in the case of Yohkoh observations dMcTiernan et al. (1999)). 
Examples of DEM analy sis with Hino de/XRT in the GOES 
C8.2 flare is presented by Ree ves et al.1 (120071) . 



2.3 Distribution of the Differential Emission Measure 



Differential emission measure, defined as DEM(T) — 



d(EM)/dT 



IdV/dT, is a useful diagnostic tool 



of non-isothermal plasma. It is usually calculated from 
a number of line fluxes formed at different temperature 
ranges. Generally the determination of the DEM distribu- 
tion is an ill-conditione d problem with no unique solution 
l&aig and Brown ( 1976b and this problem needs some kind 
of regularization. Calculations of DEM from broand-band 
filters observations is also possible if a corresponding filters 
have different temperature response functions (see Figure[T]l. 
In the case of Hinode broad-band filter observation su ch 
method was worked out and tested by I Weber et alj (|2004). 

In this paper we compare three different methods of the 
DEM reconstruction using Hinode/XRT observations with 
ten filters. The methods used are as follows: 

Genetic algorithm (GA) - This algorithm is based on 
ideas of biological evolution and natural selection mecha- 
nism. Starting from randomly chosen initial populations of 
different DEMs a new generation of DEMs is produced by 
crossover and mutations. Process of breeding (and multi- 
plication) of the whole population is controlled by assumed 
fitness criterion, e.g. DEMs with higher fitness have a 
higher probability to participate in the process of multipli- 
cation using crossover. Examples of DEM distribution cal- 
culated with genetic algor i thm c a m be found for exa mple 
in papers bvlKaastra et all dl996l) . iGuedel et all d 19971) and 



Mcintosh et al. (2000). 



We used an IDL implementation of Charbonnaues 
'PIKAIA' Fortran code Icharbonnaud d 1995b . Our popula- 
tion had 1000 individual DEM distributions or 'genes'. Af- 
ter a process of multiplication a new populations was ere- 
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Fig. 5. On the left: full Sun image obtained by XRT/Hinode with Ti__poly filter. On the right: enlarged map of NOAA 955 active region 
observed by Hinode/XKT telescope with Ti__poly filter. The DEM distribution, presented on Figure 7, was calculated for the loop pointed 
by the green arrow. 



ated consisting of 90% new DEMs and 10% of the old ones 
with the highest fitness. We used \ 2 as the fitness criterion 
of the DEM. The process of evolution was stopped after 10 
000 - 20 000 generations, when the convergence became very 
slow. The minimum values of \ 2 were around 4.3; typically 
15 to 30 'best' newly generated DEMs had \ 2 values below 
10. The evolution process was repeated 6 times, each time 
starting from a new random population. As a result we ob- 
tained 136 solutions with \ 2 < 10. In Figure|7]the green line 
presents the mean of these 136 solutions (x 2 of the mean is 
4.7). The errors were calculated as a standard deviation from 
the mean. The solution obtained with GA has \ 2 l ess than 
those obtained with WS and LS methods. 

The Withbroe - Sylwester (WS) method is iterative, mul- 
tiplicative algorithm relying on the maximum likelihood 
approach. In our calculations we started from constant 
value DEM approximation. For detailed descriptions of this 
method see Withbroe (1975) and Sylwester et al. (1980). Us- 
ing WS method we performed 200 simulations consisting of 
3000 iterations each. The obtained \ 2 were in the range from 
5 to 25. The mean solution has \ 2 = 10 (see Figure|7]). 

Least-square (LS) method - this procedure calculates 
DEM basing on a least-squares fit of the calculated fluxes 
to the observed ones. The method is described in more de- 



tail bv lGolub et all (120071) and lWeber et all d2004l) and imple- 
mented as a xrtjlemjterative routine in SolarSoft. The DEM 
function is evaluated from some spline points, which are di- 
rectly manipulated by the \ 2 &e fitting the routine. Using 
LS method we performed 200 LS-simulations. The obtained 
X 2 were in the range from 20 to 30. The mean solution has 
X 2 = 25. 

As an example we have calculated DEM for the loop ob- 



served on 13th May 2007 (see Figures [5] and |6]). The flux 
used in DEM calculation was averaged over an area delim- 
ited by 0.8 contour level on the AL_mesh filter image (see 
Figure |6l left panel). 

We have tested all three methods using synthetic isother- 
mal model. The reconstructed width of the solutions allow us 
to estimate the temperature resolutions of the method. This 
resolution is low at the edges of considered temperatures 
range (5.5 < log T < 8) and increase in the temperatures 
where response functions have their maxima (see Figure [TJ 
with best values of dlogT = 0.5-0.6, 0.2-0.3 and 0.1-0.2 for 
LS, WS and GA method respectively. DEM errors are sim- 
ilar for all three methods, and we showed in Figure [7] only 
error bars for LS for readness and transparency. 

The uncertainties are in fact very large at lower temper- 
atures (see Figure [7]), and for this reason all three methods 
give solutions lying within error limits of both other methods. 
Thus one can not evaluate precisely EM in this temperatures 
range using Hinode/XKT data. However, it is well known 
that DEM has a minimum around log T w 5.5 for quiet Sun 
and active regi ons and for this reason each obtained soluti on 
could be right dBrosius et all dl996h . lSchmelz et all d200lh ). 

There are also differences in DEM solutions obtained for 
high temperatures. The LS solution indicates the presence 
of a bigger amount of the plasma having the temperatures in 
the range (7.0 < log T < 7.3) than the WS and GA meth- 
ods. The differences are larger than estimated errors. A pos- 
sible explanation of these differences i s that LS models d e- 
pend on initial approximation of DEM (Weber et al. (12004 1. 
LS method has also the lowest temperature resolution com- 
paring to WS and GA methods. The difference between all 
three methods should be extensively tested using XRT obser- 
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Fig. 6. Images of the loop observed in active region NOAA 09555 around 18:00 UT, taken in three different Hinode/XKT filters. Contours 
overplotted on an ALjnesh filter image corresponds to 0.6 and 0.8 of the maximum flux of the loop. 
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Fig. 7. The DEM distributions of the loop observed on 13th May 
2007 obtained using the three analysed numerical methods and us- 
ing isothermal approximation with different filter pairs. 



vations. 

Obtained \ 2 are rather high. Considering whole loop in- 
stead of one or two pixels, the relative errors of flux decrease 
and thus \ 2 increase (being much above 1.0). On the other 
side there are still large uncertainties in the fluxes measured 
with Al__thick filter which strongly influence DEM in high 
temperatures. Large \ 2 value can be also attributed to uncer- 
tainness in atomic data which contribute to response func- 
tions and to systematic error in background rejection. 



3 Discussion and conclusions 

XRT telescope can detect coronal plasma of various temper- 
atures ranging from less than 1 MK up to more than 10MK. 
Applying observations made with two or more filters one can 
analyze coronal temperature and emission measure along a 
line-of-sight. 

The filter-ratio method has been used frequently in the 
past, but the high quality XRT data give temperature and 
emission measure maps with unprecedented accuracy and 
resolution. However, filter-ratio method is not appropriate 
for multi-temperature structures, giving an average value of 
T and EM only. For such a case the DEM methods of analy- 
sis of the distribution of emission measure should be used. 

The DEM distributions obtained with the described meth- 
ods are presented on Figure |7j All three methods give the 
main peak of emission for log T = 6.3. Observed differences 
between results of the three methods lie more or less within 
error limits (please note the logarithmic scale of DEM). This 
figure presents also a set of isothermal solutions obtained for 
different pairs of filters. Their values lie around the main 
peak of DEM obtained with numerical methods described 
before. Thus, we can conclude that all three methods give 
results similar to the two filters ratio method. 

We can also conclude that described methods of the DEM 
calculation gave similar and comparable solutions. Addition- 
ally, based on shapes of the emission functions (see Figure[TJ 
we can expect much better conformity of the results obtained 
with WS, LS and GA methods in a case of strong flares. We 
can state that the majority of the pairs of the Hinode filters 
allows one to derive the temperature and emission measure 
in the isothermal plasma approximation using standard diag- 
nostic based on two filters ratio method. The temperatures 
evaluated with this method for two non-flaring active region 
plasma were in the range from around 1.5 MK to about 5 
MK. 
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